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Abstract 

In this Chapter the time-domain analysis of the velocity of the electromagnetic field pulses 
generated by a spatially compact source in free space is presented. Recent simulations 
and measurements of anomalous superluminal, subluminal, and negative velocities are 
discussed. It is shown that such velocities are local and instantaneous in nature and 
do not violate either causality or special relativity. Although these effects are mainly 
confined to the near- and intermediate-field zones, some of them seem paradoxical and 
still lack adequate physical interpretation. 



1. Introduction 

Historically, the fact that the electromagnetic field propagates in vacuum with the 
speed of light c = 299792458 m/s was used to establish the electromagnetic (wave) 
nature of light and the correctness of Maxwell's equations. The invariance of the speed 
of light with respect to inertial reference frames has lead to the special relativity theory 
and its idea of space-time and was used to fix the standard of length once and for 
all, Giacomo (1984). The realization that nothing can travel faster than light was as 
profound as the deduction of the impossibility of perpetual motion from the second law 
of thermodynamics. Often, however, measurements or simulations are reported which 
seem to demonstrate velocities greater or smaller than what is expected from light. 

It has been known for a long time that the notion of velocity becomes dubious when 
a wave propagates in a medium with strong anomalous dispersion, see Landau & Lif- 
shitz (1963). The envelope of a wave in such a medium may deform during the propaga- 
tion so that it may seem to travel faster than light, slow down, completely stop, or even 
travel with a negative velocity (towards the source). Measurements confirming these 
strange effects were reported by Segard and Macke (1985), Steinberg et al (1993), Hau 
et al (1998), Wang, L. J. et al (2000), Bajcsy et al (2003), and Wang, H. et al (2007). 

With recent advances in composites, materials have been designed (notably, photonic 
crystals and metamaterials) that mimic microscopic dispersion on a macroscopic level 
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in an effective sense, i.e., their effective constitutive parameters exhibit anomalous dis- 
persion in a prescibed frequency band. There too unusual velocities were predicted and 
measured, see e.g. Spielman et al (1994) and Woodley et al (2004). 

Most surprisingly, however, strange velocities have been predicted and observed with 
free-space waves as well. For example, Bialynicki-Birula et al (2002) and Cirone et al 
(2002) discuss the evolution of the wavefunction of a free quantum particle with zero 
angular momentum. They observe that the initial propagation of the average radial 
position for a ring-shaped wave packet may be negative. This effect is specific to the 
two-dimensional case, and is related to the violation of the Huygens principle in the 
space of even dimensions. 

Recently free-space measurements of anomalous velocities with the electromagnetic 
field were reported by Ranfagni et al (1993, 1996) and Mugnai et al (2000). The authors 
have detected superluminal propagation in controlled electromagnetic beams with spe- 
cific transverse profiles. Understandably these reports steered up some discussion and 
controversy, see e.g. Hcitmann & Nimtz (1994), Diener (1996), Porras (1999), Wynne 
(2002), and Zamboni-Rached & Recami (2008). These anomalous effects, especially in 
free-space propagation, indicate an intrinsic ambiguity of the notions of velocity and 
speed when applied to the complex phenomenon of electromagnetic radiation. On the 
one hand, we have a well-defined constant of nature c = 299792458 m/s related to the 
electric and magnetic properties of vacuum as c = 1/^/eoMo- On the other hand, these 
constants show up as mere parameters in the Maxwell equations and do not directly 
refer to any kinematic process, i.e., something like a flight of a projectile. In fact, any 
kinematic analogy at all is extremely far-fetched with respect to monochromatic plane 
waves, where most of the wave-speed/velocity theory has been developed so far, and 
the notions of phase and group velocities were introduced, see e.g. Landau & Lifshitz 
(1963) for a discussion of linear dispersive media. Considerably less was done in the time 
domain, i.e. directly with pulsed signals where the kinematic interpretation would seem 
to be more natural. In this chapter we shall exploit such a time-domain approach. 

The operational notion of velocity requires a displacement of something in space dur- 
ing a certain interval of time. The ratio of the spatial displacement to the corresponding 
time interval gives the component of the velocity along a particular spatial direction, pro- 
vided that this ratio does not depend of the measurement interval. This is only true for 
the motion with a constant-velocity. Otherwise the limit of a vanishingly small interval 
must be taken, thus defining an instantaneous velocity. 

In a homogeneous medium without dispersion the light is believed to travel with a 
constant velocity. Moreover, if the source of light is localized, then we expect the wave 
to propagate strictly outwards. Hence, to measure the speed of light it is sufficient to 
place two receivers in line with the source, one further away than the other, and mark 
the time instants when the light passes each of them. Dividing the distance between the 
receivers by the difference between the detection times we should get the speed of light 
within the usual measurement errors. 

In the next section an experimental implementation of this two-point measurement 
with a short microwave pulse radiated by a small antenna in free space is considered. 
It gives a remarkable result, strikingly similar to the one reported by Ranfagni et al 
(1993, 1996) and Mugnai et al (2000). Namely, close to the source the electromagnetic 
pulse appears to travel faster than light, whereas further away the velocity of outward 
propagation approaches the speed of light. Much of the remainder of this chapter is 
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devoted to a theoretical explanation of this phenomenon from the time-domain point of 
view. We shall see that the dynamics of an electromagnetic pulse is not a simple outward 
propagation even in free space. A careful examination of the electromagnetic causality 
principle shows that this principle applies only to the relation between the source - 
current density - and the field, and to nothing else. Moreover, this relation is complex 
enough to allow not only faster than light pulses, but also pulses with apparently negative 
velocities, i.e., traveling towards the source. The latter effect was recently experimentally 
confirmed by the author, Budko (2009). 

Thus, in general, the electromagnetic field does not propagate at the speed of light. 
In fact, an instantaneous local velocity seems to be an all together more appropriate 
concept. For a pulse this local velocity approaches the speed of light from above as one 
moves further from the source. However, the wavefront, i.e., a hard to measure boundary 
between the field and the region of space which it has not reached yet, always travels at 
the speed of light. We shall also be looking into the power flow from a small transient 
source. Although, the time-integrated power flow is always positive, i.e., away from the 
source, in the near-field zone the instantaneous power may flow towards the source for 
some time. Finally, we consider the propagation velocity of the electromagnetic energy, 
defined as the ratio of the local instantaneous power flow and instantaneous energy 
density. This velocity never exceeds the speed of light in magnitude, however, may also 
become subluminal and even negative. 



2. One failed demonstration 

Outside the source, in vacuum or air, the electric field strength E(x, t) satisfies the 
homogeneous wave equation 

AE(x,t)-ls t 2 E(x,i) = Q. (1) 

c z 

This is a scalar wave equation satisfied by each Cartesian component of E individually. 
In the one-dimensional case 

d 2 x u(x,t)-^d?u(x,t) = o, (2) 

with the initial conditions 

u(x,0) = f(x), d t u(x, 0) = g(x), (3) 
the wave equation has the general (d'Alembert) solution of the form 

, \ f(x-ct) + f(x + ct) 1 f x+ct . , , 
u(x, t) = ^ '-^ >- + Y C J _ t 9( x ) dx ■ ( 4 ) 

Here one can see that for every t > the first term represents two copies of the initial 
field distribution f(x) reduced by half in magnitude and shifted in space by the amount 
±ct. This shift is proportional to c, which thus has not only the dimensions, but also the 
obvious interpretation as the speed of propagation. 
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Although things are more complicated in the three-dimensional case, the interpreta- 
tion of coefficient c as the speed of propagation follows, in particular, from the spherical- 
wave solution of Eq. ([I]) with a point-source. Hence, to demonstrate the finite and 
well-defined speed of electromagnetic waves it is sufficient to generate a current pulse 
in a relatively small antenna and show that the electromagnetic field spreads out as 
an approximately spherical wave with positive radial velocity. We use here the usual 
distinction between the velocity and the speed, where the latter is the always positive 
magnitude of the former. Throughout the text we mainly discuss a particular Carte- 
sian component of the velocity vector - along the outwards direction with respect to the 
source. In the first set of experiments this component is "chosen" by the measurement 
procedure itself. Further, when we turn our attention to the power flow, it appears to 
be the only component present in our measurement configuration (for the chosen mutual 
orientation of antennas). As we are talking about a Cartesian component of the velocity, 
its value can be positive, negative or zero. 

Modern microwave equipment allows to carry out direct time-domain two-point mea- 
surements in the classroom. For instance, the experimental setup shown in Fig. [I] was 
used in the introductory lecture of the graduate course on Electromagnetics at the Delft 
University of Technology. The author's plan was to follow the following simple routine: 

1. receive the same signal at two locations, one further away from the source than the 
other; 

2. choose some characteristic feature of the received waveform, say, the first maximum; 

3. measure the absolute arrival times of this feature at two locations; 

4. divide the difference between the distances of the two receivers to the source by 
the difference between the corresponding arrival times. 

The result of the last calculation must be equal to the speed of light in vacuum. The 
required picosecond temporal resolution is routinely achieved nowadays with an impulse 
generator producing almost identical impulses and trigger signals at a precisely controlled 
repetition rate. The remaining jitter is taken care of by averaging over many realization of 
the pulse. The absolute arrival times here are measured in the same (stationary) reference 
frame, and in fact using the same well-calibrated device - the sampling oscilloscope - 
with time-zero set at a fixed time-interval with respect to the trigger. 

Let the relative (radial) distance between the receiver locations be 12 cm. Then, the 
relative difference in arrival times must be 

f) 1 2 m 

At = nnn „' n „„ r = 0.40028 ns. (5) 

299792458 m/s v ' 

In Fig.[2]and Fig.[3]we see the results of actual measurements. The first measurement 
is performed near the source: receivers at 9 cm and 21 cm. In the second measurement 
the receivers are at 30 cm and 42 cm. Thus, both figures correspond to the 12 cm 
relative spatial separation of the receivers. As we have some freedom in choosing a 
waveform feature to mark the arrival times, we compare the results obtained using the 
first and the third maximum of the waveform. The data tips shown in the figures 
contain the absolute time values in seconds (X-axis data) at the extrema together with 
the corresponding voltage values (irrelevant to us here). The absolute time values are 
used to calculate At. 
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Figure 1: Equipment used to measure the velocity of electromagnetic waves. Upper 
photo: impulse generator on top of a sampling oscilloscope. Lower photo: dipole antennas 
attached to broadband co-axial cables (dipoles are less than 3 cm long). 

To the surprise and confusion of students only the second set of measurements, for 
the receivers at 30 cm and 42 cm, gives the correct answer, i.e., At — 0.4 ns. The 
relative difference in the arrival times closer to the source is consistently smaller, e.g. 
At = 0.37 ns in Fig. [2j showing superluminal propagation! 

Nothing can travel faster than light. Hence, either the measurements are completely 
wrong or the maxima of a waveform cannot be regarded as true physical entities, i.e., 
cannot be the carriers of information. On the other hand, such waveform features are 
routinely used to transmit data and deduce distances to objects (travel time tomography), 
see e.g. Valle et al (1999) and Skolnik (2001). Could it be that the two-point measurement 
procedure does not determine the velocity? Yet, this procedure represents the very 
definition of velocity. It turns out that the resolution of this puzzle can be found in the 
classical radiation formula relating the current in the source and the field at the receiver. 
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Figure 2: Superluminal propagation as seen in two-point near-field measurements. These 
signals were measured at 9 and 21 cm from the source (12 cm spatial separation). The 
arrival-time difference is too small: At = 0.37 ns. 



3. Causality, compatibility, continuity, and gauge fixing. 

The Maxwell equations in vacuum 

-V x H(x, t) + e a t E(x, t) = -J(x, t), 

o 

V x E(x, t) + Ai 9 t H(x, t) = 0, W 

give the local relation between the source, current density J, and the fields, electric and 
magnetic field strengths E and H. To relate the source at one location to the fields 
at another one has to solve the Maxwell equations, i.e., derive the radiation formula. 
Explicit solutions of this kind can be obtained only in some elementary cases, such as 
the homogeneous background medium we consider here. We are particularly interested 
in the way the speed of light c = enters the radiation formula, how exactly it 

reflects the causality of the radiation process, and whether this causality extends to such 
obvious and practically important features of the waveform as its extrema. 

Causality is an assumption that the fields E and H are caused by the current J. 
We must ensure therefore that the left- and right-hand sides of the Maxwell equations 
are compatible, i.e., all sources are accounted for by the current density J. To see what 
it really means let us reduce the Maxwell's system of two first-order equations to the 
second-order vector wave equation: 

V x V x E(x, t) + ^9 t 2 E(x, t) = -/i d t J(x, t). (7) 

From the first of the Maxwell equations we derive the following compatibility relation: 



e 9 t V-E(x,<) = -V- J(x,t). 
6 
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Figure 3: Further away from the source propagation happens at the speed of light. These 
signals were measured at 30 and 42 cm from the source (same 12 cm spatial separation 
as in Fig. [2j). The arrival-time difference is now as expected: At — 0.4 ns. 



Temporal causality is often introduced as an assumption that no field existed before the 
initial switch-on moment to- Hence, we can re- write the local (in time) relation ^ as a 
time-integrated formula 

1 f* 

V-E(x,t) = V- / J(x,i')di'. (9) 



For an "adiabatic" source one should take to —> —oo. Now, using this expression and 
identity V x V x E = VV • E — AE we reduce the vector wave equation Q to the usual 
scalar wave equation 



AE(x, i) — ^-<9?E(x,t) = - — VV ■ / J(x, t') dt' + /i 3 t J(x, t), 

c £o Jt 



(10) 



where the left-hand side is in its standard form (JlJ and the source term is explicitly 
identified. 

The compatibility relation ([9]), is not only a key step in transforming the Maxwell 



equations into the wave equation (10), but has a connection to gauge fixing as well. 



To see this recall that the solution of Eq. ( 10 ) can be expressed via the solution of the 



same equation with a simpler right-hand side. Indeed, transform Eq. ( 10 ) into the (k, s)- 
domain using the three-dimensional spatial Fourier and the one-dimensional temporal 
Laplace transforms: 

s 2 _ 1 - 
-(|k| 2 + ^)E(k, s) = kk • J(k, s) + s/x J(k, s). (11) 

& S£o 
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Solution of this algebraic equation can be written as 



E(k, s) = kk • A(k, s) - sA(k, s), 

s 



where 



A(M) 



/'o 



7 J(k,s), 



(12) 



(13) 



|k| 2 +S 2 /C 2 

and it is easy to deduce by inverse transformations that the new quantity A(x, t) satisfies 

1 



AA(x,i) - -d 2 A(x,i) = -MoJ(x,t). 



(14) 



The inverse transformation of Eq. (12 1 shows that there is a (spatially) local relation 



between the solutions of Eq. (10) and Eq. (14 1, namely, 



E(x,i)=c 2 VV- ( A(x,t')dt'-9iA(x,t). 

J to 



(15) 



On the other hand we know that a similar local relation exists between the fields and 
the vector and scalar potentials, i.e., 



E(x,t) = -V0(x,t) - t A(x,t). 



(16) 

We see that the two expressions, (jTHJ) and ( [Id] ), are identical, if we employ the Lorenz 
gauge 



V-A(x,t) + ia t 0(x,t)=O, 



in the following time-integrated form: 



2 V ■ f A(x,t')dt' = -0(x,t). 

Jtn 



The scalar potential, in general, also obeys a wave equation 



1 



A^(x,i)--a3ft(x,i) 



1 



p(x,t), 



which can be written as 

A^(x, t) - ^9 t 2 0(x, t) = - V • / J(x, *') d<', 
if we use the following time-integrated form of the continuity equation: 

V • f J(x,i')dt' = -p(x,t). 

Jin 



(17) 



(18) 



(19) 



(20) 



(21) 



Thus, employing the compatibility relation ([9| to derive the wave equation and the 
Lorenz gauge fixing ( 18 ) required for the unique representation of fields via potentials 



are two completely equivalent procedures. Moreover, we have also established that this 
amounts to (follows from) relating the sources of the scalar and vector potentials via the 



integrated continuity equation (21). The underlying mathematical assumption here is 



that the current density J accounts for all the sources of the electromagnetic field. Hence, 



Eq. ( 21 ) describes all possible variations of the charge density p. From the physical point 
of view, we assume that any eventual "bare" charge can only be created by some explicit 
dynamical process which disrupts the initial electrical neutrality of the substance (e.g. 
chemical reaction). Therefore, a current - motion of charges - precedes and causes the 
accumulation of charge. This is an often neglected and more subtle part of the general 
causality assumption behind the radiation formula. It treats static field as a stationary 
long-time limit of an initially transient electromagnetic field. The very existence of this 
limit and therefore static fields as such is thus an open problem in this formulation. 

4. Time-domain radiation formula 

A detailed derivation of the radiation formula presented below can be found in 
De Hoop (1995). It proceeds by computing the inverse Fourier transform of the (k, s)- 



domain solution (12 1. We arrive at 



c 2 



with 



E(x,s) = VV • — A(x,s) - sA(x,s), (22) 
s 



A(x,s) = /i / .g(x-x\ S )j(x',s)dx', (23) 
where the scalar Green's function is 

e -(s/c)|x| 

g(x,a)= , | | ■ (24) 
47r|x| 

At this stage it is convenient to carry out the spatial differentiations in the first term on 



the right in Eq. (22 1. The resulting formula is 



E ( x > s ) = I 1j I 3( Q - K] — J ( x '> s ) dx ' 

'eR 3 



e -(s/c)|x-x'| 




4tt|x-x'| 3 


Q-I] — 
se 


e -(s/c)|x-x'| 




4tt|x-x'| 2 


Q-I] — 


-(s/c)|x-x'| 

K 


8-1] A 


4tt|x-x'| 1 





+ / xi ^[3Q-I]— J(x', s )dx' (25) 

J x , eR3 4tt|x-x'| 2 e c 

+ / -n 7T [Q n J(x', s) dx'. 

Jx'GR 3 4tt|x-x'| £qC Z 

The three terms above are called near-, intermediate-, and far-field contributions in 
accordance with their spatial decay factors. The tensors act as 



x — X' 

J(x')=J(x'). (27) 
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The magnetic field can be found from a similar formula 



r P -( s / c )l x - x I 
H(x, S ) = -/ — — 0x J(x', S )dx' 

e -(«/c)|x-x'| 1 * 

— -0 x J(x',s)dx', 



x'ei 



where 



47r|x — x'l c 



Finally, the inverse Laplace transform gives the explicit time-domain formulas 

E(x, t)= [ - i^' 1 ,,.. - J(x', t 1 ) dt' dx' 
7 x , eR 3 47r|x-x'|3 £o J t0 

+ / , , Q ~\ 3 — J(x',^,)dx' (30) 

Jx'GR3 47T|X-X'|^ £ C 



and 



H(x -"^X*J^T exJ(lt ' ili,dx ' 

-( ,tebr^i ex ? 3 " ,(x '- <B) ' dx ' 

ix'fB 3 47r l x ~ x I c 



(31) 



/x'G 

where for the first time we meet the retarded time 



tR = t- ] ±^±, (32) 



which is an analogue of the ±ci shift of the d'Alembert solution Q. The spatial inte- 
gration in (30 1 will in practice be limited to a finite domain D s occupied by the source. 



5. Explanation of near-field superluminal velocities 



With the radiation formula ( 30 1 at hand we can start analyzing the shape of the 
waveforms received at some distance from the source. A waveform measured by the 
dipole antenna of the type shown in Fig. [T] gives the time-evolution of a single component 
of the electric field strength determined by the orientation of the dipole and modified by 
the the antenna-cable-oscilloscope receiving tract. The latter modification, however, is 
the same for all receiver locations, i.e., does not depend on the distance from the source. 
Only in the close proximity of the source the expected small multiple reflections between 
the two antennas may slightly modify the tail of the waveform with respect to what is 



predicted by Eq. ( 30 1 
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Figure 4: Back-in-time shifting of extrema. Depicted are the signal, its cumulative time- 
integral, and derivative. The first extremum of the integral is further (later) in time, 
whereas the first extremum of the derivative is earlier in time than the first extremum 
of the original signal. 

A simplified formula for small source and receiver dipoles, parallel to each other and 
located in a plane orthogonal to their orientation (as in Fig. [IJ is 

V ^ *> ~ "4^ C m d ' ~ J^ I{tR) ~ ^R dtI ^ (33) 

where R is the distance between the dipoles, in . = t — R/c, V is the measured waveform, 
and / is the signal (current) in the source antenna. 

Obviously, the received waveform is a function not only of the signal I at an earlier 
time, but also of its time-derivative and time-integral. Figure [4] shows the three functions 
in question for the current of the form 

I(t) = -2a(i-t c )fr Q ( 4 -^ 2 , (34) 

whose integral and derivative can be computed analytically. As can be seen from Fig. [4] 
the extrema of the three functions are time-shifted with respect to each other. In the 
present case we can talk about the first extremum only, since the integral does not have 
a second one. With respect to this first extremum we notice that the integral has it 
later, and the derivative has it earlier in time than the function itself. This is a purely 
mathematical phenomenon and it will be observed with any function. 

The received waveform is a weighted sum of these three functions. The weight coef- 
ficients are the functions of distance R between the source and the receiver, where the 
term with R~ 3 will be dominant close to the source (near- field zone), and the term with 
is dominant for large R (far-field zone). This means that in the near-field zone the 
received waveform will look almost like the integral, and in the far-field zone it will look 
almost like the derivative of the current. 
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The two-point procedure described in Section 2 deduces the velocity of the electro- 
magnetic pulse from the waveform measured at two locations where one is further away 



from the source. Thus, according to Eq. (34) the values of the weighting factors will be 



different for these two locations. Not only the waveform will be smaller in amplitude for 
the more distant of the two locations, but the relative weight of the three contributions 
will change as well. In general the relative weight of the last term increases, whereas the 
relative weight of the first term decreases with distance. Hence, we may expect the shape 
of the waveform to gradually change from the one dominated by the time-integrated sig- 
nal to the one dominated by the time-derivative of the signal. As can be seen from Fig. [4] 
this means a shift of an extremum towards earlier times. Such backward in time "mo- 
tion" will be superimposed on the normal time-delay which shifts the received waveform 
to the right along the time axis as one moves further away from the source. 

The latter observation explains our failure to measure the speed of light in Section 2. 
Indeed, the first measurement location in Fig. [2] is close enough to the source for the 
influence of the near- and intermediate-field terms to be significant. Therefore, the 
overall waveform is somewhat shifted to the right along the time-axis relative to the 
time-delayed derivative of the original signal, which is dominant at larger distances. At 
the second location the influence of the near- and intermediate-field terms is smaller and 
so is the relative shift. Let us express the arrival times of the measured extrema at two 
locations as 



R 

tarr(-Rl) = *0 + — + Af(fli), 
C 

fan(J2a) = *o + — + Ai(-R 2 ), 
c 



(35) 



where At(R) is the relative shift due to the near- and intermediate-field terms. We give 
it a positive sign to emphasize that the relative time-shift is happening towards later 
(positive) times on top of the normal propagation induced positive time-delay R/c. 

The speed of light was deduced in Section 2 from the following simple two-point 
calculation: 

R2 ~ Rl (36) 



tarr{R2) ~ t gxt (Rl) ' 

applied to the extrema of the measured waveform. It is easy to see, that without the 



additional relative time-shifts At, substitution of (35) would have given us the exact 



speed of light. Now, however, the denominator in (36) is 



tarr(#2) " t«r(#l) = ~ ^ + Mifc) - At^). (37) 

C 

Since the relative time-shift At(R) is diminishing with distance, we have in general 

At(R 2 ) - At(i?i) < 0, R 2 > Ri, (38) 

Thus, the measured time difference is always smaller then the expected one. If this time 
difference is positive we get 



Ccxp C, 
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(39) 



i.e., the measured speed is superluminal - exactly as we have observed. 

That this superluminal behavior is not visible in Fig. [3j and in the majority of exper- 
iments with light, can also be explained. The near- and intermediate-field terms, which 
cause the additional time shift, decay not only relatively with respect to the far-field 
term, but also absolutely, i.e., their influence is practically not measurable in the far 
field. The far-field zone is a notion relative to the wavelength. Usually, the far-field zone 
starts beyond a few tens of wavelengths, i.e., is very close to the source for the visual 
light. 

The measured speed reflects two competing phenomena - outward propagation and 
the diminishing additional positive time shift. In the far-field zone the first of these 
phenomena gradually takes over, and we may expect the local two-point measurement 
procedure to yield values of speed progressively approaching the speed of light from 
above, as if the pulse was decelerating. In the near-field zone, however, the additional 
time shift is large enough to be detectable. We saw how it yields superluminal values 
for the pulse speed. Obviously, we need to let go of the idea of a constant speed and 
consider a more general concept of a local velocity which is a function of distance from 
the source. For example, we could introduce a local velocity as a limit of the two-point 
measurement procedure in the neighborhood of some location R as 

vioc(R) = hm — — ^ — — , (40) 

where t ex (R) is the time of the waveform extremum measured at location R. This 
definition, however, is hard to apply as we do not know the above extrema times explicitly. 



6. Locally negative velocities 

As was mentioned above superluminal pulse propagation is observed, if the measured 
time difference in the two-point procedure is smaller than the expected R/c, but remains 
positive. If this difference is negative, then the measured velocity would also become 
negative. Are such negative velocities possible, and what does it mean? The answer 
can be obtained with the formulas and understanding we have acquired in the previous 
section. Indeed, a negative velocity could be measured in a two-point procedure, if 
the anomalous positive time-shift we discussed was disappearing faster with R than the 
normal growth of the positive time-shift due to propagation, i.e., 

R * ~ Rl < At(Rt) - At(R 2 ). (41) 
c 

Recently an experiment was reported by the author, Budko (2009). It showed that this 
effect is indeed present and can be observed with the setup which we used above to 
demonstrate the superluminal pulse velocities. Very close to the source the waveform 
shows a small but detectable movement towards earlier times (to the left on the scope 
screen) as we move the receiver away from the source in small steps. Figure [5] shows the 
results of this experiment. The measurements had to be carried out very close to the 
source. For the 75-ps pulse that was used in the experiment the negative velocity region 
is within 10 mm distance from the source. The extent of the negative- velocity region was 
initially predicted by a simulation similar to the one considered in the next section. The 
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Figure 5: Experimental demonstration of local negative pulse velocity in the near-field 
of a small dipole antenna. Upper image shows the space-time dynamics of the signal. 
Close to the source (bottom of the image) the wave makes a small bend to the left along 
the time axis - negative velocity region. Further from the source the wave dynamics 
approaches the light-cone behavior (white line). Lower plot shows the wave dynamics 
in the near-field zone. While the wavefront and the first slope shift towards later times 
(positive velocity), the extrema and the whole inner part of the waveform shift to the 
left (negative velocity). 

actual measurements also contain multiple reflections between the source and receiver 
dipoles. These reflections show as a small bump in the slope of the waveform between 
cot = 0.1 and cgt = 0.15 m. This was verified by inserting a larger metallic scatterer 
in the plain of the source and observing an increase of the said bump as well as the 
appearance of multiples further down the tail of the waveform. Apart from that the 
leftward shift of the inner part of the waveform is clearly visible. The upper image in 
Fig. [5] is a rotated version of a typical light-cone diagram from special relativity. In fact, 
the white line indicates the way the light cone is supposed to be. If the light was indeed 
following its "cone", then the colored stripes in the image should all follow that white 
line. This happens only further away from the source. Whereas in the bottom of the 
image there is an obvious near-field distortion - a leftward bend corresponding to the 
negative velocity of propagation. 

Such behavior is very counterintuitive and may even seem anti-causal. Indeed, it 
looks like the signal reaches a more distant observer faster than a closer one. Although 
information may be and often is encoded in extrema of a waveform the speed at which it 
reaches the observer is hard to define. We notice from the experimental data in Figure [5] 
that the front of the wave travels with a positive velocity. Moreover, for a source current 
with a sharp temporal boundary I(t) = 0, t < t we can prove analytically that the 
wavefront always travels outwards at the speed of light, since near- and intermediate- 
field terms are exactly zero for t — £r in that case. For the same reason the end of 
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the waveform would also travel outwards and at the speed of light. Thus, the beginning 
and the end of an information carrying wave-packet would arrive at two observes in a 
proper relativistically consistent temporal order, i.e., later for the more distant observer. 
However, due to the anomalous temporal shift of extrema the more distant observer may 
actually receive the information sooner than the closer one. Of course, we are talking 
about the near-field zone and two very close observers, but nevertheless. Explanation 
seems to be in the fact that the two observes are causally related to the source of the 
field, but not to each other. The situation becomes less paradoxical if we put it this 
way: it is not that the more distant observer gets the information faster than normal. 
No, the beginning of the information packet arrives at the speed of light. It is that 
the closer of the two observers receives the same information slower than normal, due 
to the additional positive time-shift of the waveform in the near-field zone. The same 
explanation helps to comprehend the superluminal impulse velocities observed earlier. 



7. Negative power flow 



At the moment the strange velocities we observe in the near-field zone seem to be 
tightly linked to the two-point measurement procedure and our choice of extrema as 
the reference points. To show that the two-point procedure is not really essential here 
let us look at the power flow in the near-field zone. An analysis along the following 
lines, although focused on the Hertzian dipole with an initial static field, can be found 
in Schantz (2001). The instanteneous power density flowing through location x at time 
t is given by the Poynting vector 



S(x,t) = E(x,t) x H(x,t). 



(42) 



In our particular experimental arrangement with a small source the current density can 
be approximated by 



J(x,t) = d/(i)<5(x-x s ), 



(43) 



where d is a unit vector in the direction of the current. For the geometry depicted in 
Fig. [T] we have 



S(R,t) = P(R,t)&, 



(44) 



where the unit vector points away from the source towards the receiver and the power 
density amplitude is 
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Figure 6: Top: mono-cycle source current, its integral, and derivative (normalized). 
Bottom: Simulated electric field at various distances from the source (normalized) . Close 
to the source the received signal resembles the integral taken with a negative sign. Further 
it deforms into a shifted negative of the signal, and finally transforms into a negative of 
the derivative of the signal. The negative velocity phenomenon is obvious between 1 mm 
and 15 mm distances - waveform shifts leftwards. 



While the last term is undoubtedly positive at all times, the remaining two terms may, 
in principle, become negative due to the mixed products of the current, its derivative, 
and integral, which as we know can all have different signs at the same time - see Fig. [4] 
If a negative power flow can be detected, the phenomenon would again be confined to 
the near-field zone, as the dominant far-field term is positive. 

Let us consider another current pulse with analytically known derivative and integral. 
This time, however, we shall look at a pulse called monocycle with well-defined beginning 
and end: 

TM - I sin ( wt )' te[t ,t + 2n/w], . . 

0, t^[t ,t + 2n/uj] ) l4bj 

where to is the carrier frequency. The derivative and integral of this signal within the 
time interval [to, to + 2tt/uj] are: 

d t I(t) — wcos(wi), 

(47) 



1(1/) dt' = - [cos(wto) - cos(w*)] , 

and are zero outside the said interval. The normalized versions of these functions are 
shown in the top plot of Fig. [6] for ui = 8n x 10 9 , i.e., for a 4 GHz carrier frequency. The 
bottom plot of that figure shows the results of simulation of the electric field at three 



different distances from the source: 1, 15, and 50 mm. We use Eq. (33) to simulate the 
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Figure 7: Normalized time-integrated power as a function of distance from the source. 
This quantity is positive everywhere. 



waveforms. All pertaining deformations of the waveform are now clearly visible, as it 
resembles the integral of the signal in the near-field zone and gradually transforms into 
the derivative of the signal as we measure further away. We also see the negative velocity 
phenomenon as the inner part of the waveform at 15 mm is clearly shifted to the left, 
instead of being shifted to the right of the waveform corresponding to 1 mm distance. 
At the same time the wavefront and the end of the waveform do shift rightwards as 
predicted. 

To understand, at least mathematically, the emergence of the negative power flow 
consider a very small distance R, such that the first terms in Eq. (45) are dominant. 
These terms contain mixed products of the current, its integral, and derivative. For 
the particular current given by Eq. (46)— (47) we see that these mixed products will be 
negative over a certain period of the time - see Fig. [8] (top). Although we do not prove it 
here rigorously for an arbitrary current, computer simulations and common sense indicate 
that the time-integrated power is positive at all distances, including the near-field zone - 
see Fig. [7j This, however, does not exclude the possibility of instantaneous power being 
negative for some time. 

In Figure [8] (bottom) we see the power P(R : t) as a function of time at three distances 
from the source. The negative power flow is clearly visible in this plot. As we move 
towards the far-field zone the power flow looks positive at all times. This, however, is 
not true as we shall see in the next Section. 

Power flow tells us how much energy is coming from a certain direction in a unit of 
time. If a power flow is negative it means that the energy is decreasing with time, or that 
the power is flowing towards the source of radiation, instead of away from it. Since the 
Poynting vector is related to the density of the linear momentum of the electromagnetic 
field, we may expect the instantaneous pressure of light to change its direction over 
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Figure 8: Top: normalized products of the signal and its integral, and derivative of 
the signal and signal's integral corresponding to the source current shown in Fig. [6] 
Notice the presence of negative parts in both cases. Bottom: normalized instantaneous 
power flow at various distances from the source. In the near-field zone the power flow 
can become negative for some time. During that time power flows towards the source. 
Further from the source the power flow looks positive at all times (see, however, Fig.[9|. 



extremely short time intervals in the near-field zone. Not only the light will be pushing 
away from the source, as it is usually the case, but pulling towards the source as well. 
This may induce a very rapid mechanical oscillation on a neutral test particle - an opto- 
mechanical effect which might be of interest in nano-physics. 



8. Energy velocity 

We can formally relate the power flow and the energy density in our experimental 
setup via the following formula: 

P(x, f)0 = v (x, t)£ (x, t)0, (48) 

where £ is the local instantaneous energy of the electromagnetic field. From where the 
energy velocity along is expressed as 

Ua (x,t) = £M. (49) 
£{x,t) 

Simulations of the energy velocity for a transient Hertzian dipole were presented by 
Marcano & Diaz (2006). Since, as we saw above, the local power flow can be negative, 
while the field energy is always positive, the near-field zone results of Marcano & Diaz 
(2006) naturally contained locally negative energy velocities. 
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Figure 9: Top: normalized local instantaneous energy at various distances corresponding 
to the source current shown in Fig. [6j Bottom: Relative energy velocity, v c /c, at various 
distances from the source. 



The power amplitude for our measurement configuration is given by Eq. (45 ), whereas 
the energy density at location R as a function of time is 
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This function for the monocycle source current Eq. ( 46 1 and the corresponding instan- 
taneous relative energy velocity v e (R,t)/c are shown in Fig. [9} Similar to the waveform 
of the electric field, the waveform of the instantaneous energy deforms as one moves 
through the near-, intermediate-, and far-field zones. There is no doubt about the posi- 
tivity of energy density at all locations and times due to the intrinsic mathematical form 
of expression Eq. (50 1. Thus, the negative values of the energy velocity visible in Fig. [9] 
(bottom) are caused by the locally negative power amplitude. It is obvious, that the 
energy velocity varies between — c and c. It is also obvious that the energy velocities of 
the front and the end of the wave-packet are both exactly the speed of light, and that the 
deviations concern the inner part of the waveform only. We notice that the subluminal 
and negative velocities are mainly featured in the near-field zone and that the far-field 
energy waveform moves mostly at the speed of light. However, at the time instant where 
both energy and power tend to zero, the energy velocity makes a sharp swing to — c and 
back even in the far-field zone. This shows that the power amplitude must be going 
slightly negative at the minima of the far-field waveform, something we could not previ- 
ously detect in Fig. [HJ Thus, although the local instantaneous energy velocity seems to 
be a very natural concept, and its magnitude is explicitly bounded by the speed of light, 
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it also magnifies the effect of the first two terms in Eq. ( 45 1 , so that it is felt even in the 
far-field zone, albeit for a vanishingly small time. Whether this swing in the power flow 
can be detected remains to be seen. 



9. Conclusions 

We have discussed a series of unusual transient effects in the electromagnetic radiation 
from a localized source. These effects demonstrate that the space-time evolution of 
the electromagnetic field is nothing like a simple outwards propagation of a spherical 
wave at the speed of light. In fact, the speed of light could only be assigned to the 
boundary between the field and the area free from it - the wavefront and the end of 
a wave-packet. The inner part of the wave-packet undergoes significant deformation in 
the course of propagation. This deformation leads, in particular, to superluminal and 
negative velocity results with the classical two-point velocity measurement procedure. We 
have uncovered that the source of this deformation is an anomalous time delay, which 
the waveform undergoes in the immediate vicinity of the source. It takes more time 
for, say, an extremum to build up close to the source than a little further away. Hence, 
paradoxically, it takes more time for the information to arrive at a closer observer than 
at an observer further away from the source. Yet, both transmission times are within 
what is allowed by the relativity theory. While this effect can be observed only very close 
to the source, the other related phenomenon can be measured in the intermediate-field 
zone as well. 

The anomalous time delay rapidly decays with the distance from the source. Hence, 
the time-interval measured in the two-point procedure is usually smaller than what is 
expected from the propagation at the speed of light. This explains the observed su- 
perluminal results and shows the necessity of near-field corrections to the travel-time 
algorithms used in geophysics and radar imaging, such as those used in e.g. Valle et al 
(1999) and Skolnik (2001). 

We have also considered the instantaneous power flow from a transient source. In 
the near-field zone power may flow towards the source for some time. Although, the 
time-integrated power flow stays positive, a positive instantaneous power flow followed 
by a negative flow during a significant interval of time will induce an oscillating linear 
momentum and an oscillating instantaneous light pressure. 

Finally, we have analyzed the local instantaneous energy velocity which is conve- 
niently bounded between — c and c. Its waveforms show that the energy velocities of the 
wavefront and the end of a wave-packet are both exactly c, and that the subluminal and 
negative values during extended time intervals are the feature of the near-field zone only. 
At the same time the energy velocity reveals that the power flow can become negative 
for a short time even in the far-field zone. This happens when both the power flow and 
the energy are close to zero, and it remains to be seen whether this effect can be detected 
exp eriment ally. 

On a theoretical level all these unusual effects and their experimental confirmation 



show the correctness of the time-domain radiation formula (30) -(31). We have demon- 
strated how this formula follows from the Lorenz gauge-fixing between the classical scalar 
and vector potentials and eventually reduces to the continuity equation between the cur- 
rent and charge densities. These otherwise somewhat ad hoc procedures thus acquire 
additional physical significance. 
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